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Abstract. We give a proof of stability of relativistic matter with magnetic fields 
\Q ' all the way up to the critical value of the nuclear charge Za = 2 /it. 
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. 1. Introduction 

o, 

, We shall give a proof of the 'stability of relativistic matter' that goes further than 

previous proofs by permitting the inclusion of magnetic fields for values of the nuclear 
charge Z all the way up to Za = 2/ir, which is the well known critical value in 
the absence of a field. (The dimensionless number a = e 2 /he is the 'fine-structure 
\q ; constant' and equals 1/137.036 in nature.) More precisely, we shall show how to 

modify the earlier proof of Theorem 2 in |LY| so that an arbitrary magnetic field can 
be included. Reference will freely be made to items in the |LYj paper. 

o 

of stability of matter, will be given in the next section. For a detailed overview of this 
topic, we refer to |L1[ IL2j . For the present we note that stability requires a bound on 



a in two ways. One is the requirement, for any number of electrons, that Za <2/n. 
In fact, if Za > 2/n the Hamiltonian is not bounded below even for a single electron. 
The other requirement is a bound on a itself, a < a c , even for arbitrarily small Z > 0, 
which comes into play when the number of particles is sufficiently large. It is known 
that a c < 128/1 5tt; see jlYl Thm. 3] and also jDLj . 

For values of Za strictly smaller than the critical value 2/7T, it has been shown that 
stability holds with a magnetic field included. This is the content of Theorem 1 in 
|LYj . in which the critical value of a c goes to zero as Za approaches 2/tc, however. 
(The result in |LY| Theorem 1] does not explicitly include a magnetic field, but the 
fact that the proof can easily be modified was noted in |LLoSoj .) A similar result, by 
a different method, was proved in |LLoSij . 

The more refined Theorem 2 in |LY| gives stability for the 'natural' value Za < 2/tt 
and all a < 1/94. While the true value of a c is probably closer to 1, the value 
1/94 > 1/137 is sufficient for physics. The problem with the proof of |LY[ Theorem 
2] is that it does not allow for the inclusion of magnetic fields. Specifically, Theorems 
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9-11 have to be substantially modified, and doing so was an open problem for many 
years. This will be accomplished here at the price of decreasing a c from 1/94 to 1/133. 
Fortunately, this is still larger than the physical value 1/137 ! 

In a closely related paper [FLSeJ we also show how to achieve a proof of stability for 
all Za <2/tt with an arbitrary magnetic field, but the value of a c there is very much 
smaller than the value obtained here. In particular, the physical value of a = 1/137 is 
not covered by the result in |FLSej . The focus of |FLSej is much broader than 'stability 



of matter', however. It is concerned with a general connection between Sobolev and 
Lieb-Thirring type inequalities, and includes as a special case Theorem 14.51 of this 
paper. The proof of the general result in jFLSej is much more involved than the one 
of the special case presented here, and yields a worse bound on the relevant constant. 

Acknowledgements. We thank Heinz Siedentop for helpful remarks. This work 
was partially supported by the Swedish Foundation for International Cooperation in 
Research and Higher Education (STINT) (R.F.), by U.S. National Science Foundation 
grants PHY 01 39984 (E.L.) and PHY 03 53181 (R.S.), and by an A.P Sloan Fellowship 
(R.S.). 

2. Definitions and Main Theorem 

We consider N electrons of mass m > with q spin states (q = 2 for real elec- 
trons) and K fixed nuclei with (distinct) coordinates Ri, . . . ,Rk € M 3 and charges 
Zi, . . . , Zr > 0. The electrons interact with an external, spatially dependent magnetic 
field B{x), which is given in terms of the magnetic vector potential A(x) by B = curlA 
A pseudo-relativistic description of the corresponding quantum-mechanical system is 
given by the Hamiltonian 

H n ,k ■= ^2 ( \J (Pi + A ( x j)) 2 + m 2 - mj + aV N}K (x 1 , . . . ,x N ;Ri, . . . ,R k ) ■ (2.1) 

The Pauli exclusion principle for fermions dictates that H^ k acts on functions in the 
anti-symmetric iV-fold tensor product A N L 2 (R 3 ; C q ). We use units in which ft = c = 1, 
a > is the fine structure constant, and 

N K 

Vn,k( x i, ■ ■ ■ j %n] Ri-, • • • Rk) '•— \xi — Xj\ 1 — Z k \xj — Rk\ 1 

l<i<j<N 3=1 k=l 

+ Yl Z^Rk-Ri]- 1 . (2.2) 

l<k<l<K 

is the Coulomb potential (electron-electron, electron- nuclei, nuclei-nuclei, respectively). 
In this model there is no interaction of the electron spin with the magnetic field. Note 
that we absorb the electron charge y/a into the vector potential A, i.e., we write A(x) 
instead of \faA(x) in (|2.ip . Since A is arbitrary and our bounds are independent of 
A, this does not affect our results. 
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Stability of matter means that Hn,k is bounded from below by a constant times 
(N + K), independently of the positions Rk of the nuclei and of A. For a thorough 
discussion see |L1| IL2j . By scaling all spatial coordinates it is easy to see that either 
infij fc) ^(inf spec Hn,k) = —mN or = — oo. 

We shall prove the following. 

Theorem 2.1 (Stability of relativistic matter with magnetic fields). For qa < 

1/66.5 and aZj < 2/tt for all j , 

Hn,k > —mN 

for all N, K, Ri, . . . , Rk and A. 

For electrons q = 2 and hence our proof works up to 

1 1 

a = > . 

133 137 

The rest of this paper contains the proof of Theorem 12.11 but let us first state an 
obvious fact. 

Corollary 2.2. As a multiplication operator on A N L 2 {R 3 ; C q ), 

N 

Vn,k(%i, ■ ■ ■ ,%n)Ri, ■ ■ -Rk) > — max{66.5g, ttZj/2} } j \pj + A(xj)\ (2.3) 
for all A. 

This, of course, is just a rewording of Theorem 12. 1[ but the point is that it provides 
a lower bound for the Coulomb potential of interacting particles in terms of a one- 
body operator \p + A(x)\. This operator is dominated by the nonrelativistic operator 
\p + v4(x)| 2 and, therefore, ()2.3|) is useful in certain nonrelativistic problems. For 
example, an inequality of this type was used in |LLoSoj to prove stability of matter 
with the Pauli operator \p + A(x) | 2 + a ■ B(x) in place of \p + A(x) | 2 . It was also used 
in jLSiSoj to control the no-pair Brown- Ravenhall relativistic model. 

An examination of the proof of Theorem 2 in |LY| shows that there are two places 
that do not permit the inclusion of a magnetic vector potential A. These are Theorem 9 
(Localization of kinetic energy - general form) and Theorem 11 (Lower bound to 
the short-range energy in a ball). Our Theorem 13.11 is precisely the extension of 
Theorem 9 to the magnetic case. It may be regarded as a diamagnetic inequality on 
the localization error. It implies that Theorem 10 in |LY| holds also in the magnetic 
case, without change except for replacing \p\ by \p + A\; see Theorem 13.21 below. 

A substitute for Theorem 11 in |LYj will be given in Theorem 14.51 below. It is based 
on the observation that an estimate on eigenvalue sums of a non-magnetic operator 
with discrete spectrum implies a similar estimate (with a modified constant) for the 
corresponding magnetic operator. This is not completely obvious, since there is no 
diamagnetic inequality for sums of eigenvalues. (In fact, a conjectured diamagnetic 
inequality actually fails for fermions on a lattice and leads to the 'flux phase' |L3j.) It 
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is for the different constants in Theorem 11 in |LY| and in our Theorem 14.51 that our 
bound on a c become worse than the one in |LYj . 

As should be clear from the above discussion, our main tool will be a diamagnetic 
inequality for single functions. The one we use is the diamagnetic inequality for 
the heat kernel. In the relativistic case it states that for any A G L^ oc (IR 3 ; M 3 ) and 
u G L 2 (M 3 ) one has 

\(exp(-t\p + A\)u)(x)\<(exp(-t\p\)\u\)(x), x G R 3 . (2.4) 

This follows with the help of the subordination formula 

151 _ f°° e -t-|€l7(«) dt 



Tit 



from the 'usual' (nonrelativistic) diamagnetic inequality for the semigroup exp(— t\p- 



A\ ); see, e.g., |S3j . The heat kernel is not prominent in [LYJ, and our reformulation 
of some of the key estimates in |LY| in terms of the heat kernel is the principal novel 
feature of this paper. 

3. Localization of the kinetic energy with magnetic fields 

3.1. Relativistic IMS formula. In this subsection we establish the analogue of The- 
orem 9 in |LYj in the general case A ^ 0. First, recall that the IMS formula in the 
nonrelativistic case says that for any u and A 

/ \(p + A)u\ 2 dx = V / \( p + A)( X] u)\ 2 dx- / V|V Xj |V| 2 dx 

whenever Xj are rea l functions with Y^=iXj = 1- In this case the localization error 
Sj=o I^X? I 2 ^ s l oca l an d independent of A. The analogue in the relativistic case is the 
following special case of |FLSet Lemma B.l]. For the sake of completeness, we include 
its proof here. 

Theorem 3.1 (Localization of kinetic energy — general form). Let A G L^ oc (IR 3 ; M 3 ) 
// Xo, ■ ■ ■ , Xn are real Lipschitz continuous functions on IR 3 satisfying Y^j=o Xj = 1; 
then one has 

n 

(u, \p + A\u) = ^{XjU, \p + A\xju) - (u, L A u) . (3.1) 

j=0 

Here La is a bounded operator with integral kernel 

■n 

L A (x, y) := k A (x, y) ^(Xj(«) - Xj{y)f , 

j=0 

where k A (x,y) := lim t | t~ l exp(— t|p + v4|)(x, y) for a. e. x,y G M 3 and 
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Note that ()3.2|) says that 

1 - 

\L A (x,y)\ < L(x,y) := ~ ^j(^)) 2 • ( 3 - 3 ) 

Z7r l x f I i=0 

Here, L(x,y) is the same as in |LYt Eq. (3.7)]. Therefore, (|3.2j) is a diamagnetic 
inequality for the localization error. 

Proof. We write k A (x,y,t) := exp(— t\p + A\)(x,y) for the heat kernel and find 



J2iXjU, (1 - exp(-t|p + A\)) Xj u) = (u, (1 - exp(-t\p + A\))u) 
\i>2jj k A(x,y,t)(xj(x) - Xj(y)fu{x)u{y) dxdy. 



3=0 

+ 



(This is proved simply by writing out both sides in terms of k A (x,y,t) and using 
^2 X 2 j = 1-) Now we divide by t and let t — > 0. The left side converges to Y^j=o{Xj u i \p+ 
A\xju). Similarly, the first term on the right side divided by t converges to (u, \p+A\u). 
Hence the last term divided by t converges to some limit (u,Lau). The diamagnetic 
inequality (|2.4j) says that 

\k A {x,y,t)\ < exp(-t\p\)(x,y) = —— — 2 

tt z (\x — y\ z + t z ) 

(see |LLo| Eq. 7.11(9)]). This implies, in particular, that La is a bounded operator. 
Now it is easy to check that La is an integral operator and that the absolute value of 
its kernel is bounded pointwise by the one of L in (j3.3j) . ■ 

3.2. Localization of the kinetic energy. In this subsection we will bound the 
localization error La by a potential energy correction and an additive constant. This 
is the extension of Theorem 10 in jLYj to the case A ^ 0. It is important that both 
error terms in our bound can be chosen independently of A. 
First we need to introduce some notation. We write 

Br := {x : \x\ < R} 

for the ball of radius R and xb r f° r its characteristic function. If R = 1, we omit the 
index in the notation. We fix a constant < a < 1 and Lipschitz continuous functions 
Xo, Xi with xl + Xi = 1 such that suppxi C B\- a . With these we define L as in ()3.3j) 
with n = l. We decompose L in a short-range part L° and a long-range part L 1 given 
by the kernels 

L l (x, y) := L(x, y)xi3(x)xB(y)XB tJ ( x ~ y), L°(x, y) := L(x, y) - L l (x, y). (3.4) 
Define 

n:=lTr(L ) 2 (3.5) 
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and, for an arbitrary positive function h on B, 

9{x) := h' x {x) / L 1 (x,y)h(y)dy = h~ 1 (x)xB(x) / L(x,y)h(y) dy . 

JB J\y\<l,\x—y\<a 

Finally, for e > we define the function 

U* £ :=e X B^ + e (3.6) 
and note that U* is supported in B. 

Theorem 3.2 (Localization of kinetic energy — explicit bound in the one-center 
case). For any e > and any non-negative trace-class operator 7 one has 

1 

Tr 7 bp + A\ > J2^XjlxA\P + A \~ U;) - e-^hll- (3.7) 

j=0 

For A = this is exactly Theorem 10 in |LYj . As explained there, U* is a potential 
energy correction with only slightly larger support than Xi- The last term in (|3.7J) is 
due to the long range nature of |p + A|. It depends on 7 through its norm ||7|| but not 
through its trace. We emphasize again that both error terms in the inequality (|3.7J) 
are independent of A. 

Proof. The localization formula (|3.1|) yields 

1 

Tr7|p + A\ = XjlXjlP + A\ — Tr 7^, 

3=0 

so we only have to find an upper bound for Tr7L^. We decompose La = L° A + L\ in 
the manner of ()3.4j) and, following the proof of Theorem 10 in |LYj word by word, we 
obtain 

Tr 7 L° < eTr 7 X Bl „ CT + (2£)- 1 || 7 || Tr (L°) 2 , Tr 7 ^ < Tr 7 ^- 
Here 9a(x) := if x ^ B and, if x G B, 

e A (x) := h^(x) / \L A (x,y)\h(y)dy. 
Jb 

The estimate \La(x, y)\ < L(x,y) from Theorem 13.11 implies that Tr (L A ) 2 < 2Vt and 
that 9a < 9. This leads to the stated lower bound. ■ 

4. Bounds on eigenvalues in balls 

So far we have considered |p + A\ and its heat kernel. Now we address \p + 

A\ — 2/(7r|x|) and its heat kernel. First of all, let us recall Kato's inequality [K&; 
Eq. (V.5.33)] 

(u,\p\u)> {2/tz){u,\x\- 1 u). (4.1) 

(See also piWI lKPSj .) 
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Now let r C R 3 be an open set (we shall be interested in the case where T is a ball) 
and consider the quadratic form given by Qr(u) = (u, (\p\ — 2/tt\x\)u), restricted to 
those functions u G L 2 (M 3 ) that satisfy u = on T c , the complement of T. Of course, 
we also require u to be in the quadratic form domain of \p\ — 2/ir\x\. The quadratic 
form Qr is non-negative by (j4.1)l and it is closed (because the form \p\ — 2/tt\x\ is 
closed on L 2 (R 3 ) and limits of functions that are zero on T c are zero on T c ). From this 
it follows that there is a self-adjoint operator H-p on some domain in L 2 (T) such that 
Qr(u) = (u,Hpu). With this operator, we can define the 'heat kernel' exp (—tH r ) on 
L 2 (r) and its trace. (The fact that the trace is finite when the volume of T is finite 
follows from subsequent considerations.) 

Similarly, for a magnetic vector potential A 6 L 2 oc (R 3 ; IR 3 ), we define the operator 
in L 2 (T) using the quadratic form (u, (\p + A\ — 2/n\x\)u). Note that ()2.4|) implies 

that 

(u,\p + A\u) > {\u\,\p\\u\) . (4.2) 
This, together with ()4.1|) . shows that (u, (\p + A\ — 2/ir\x\)u) is non-negative. 

Lemma 4.1 (Heat kernel diamagnetic inequality). Let T C M 3 and let A e 

L 2 oc (M 3 ; E 3 ). Then, for any t > 0, 

Tr L2(r) exp (-tH£) < Tr L2(r) exp (-tH r ) . (4.3) 

Proof. For n = 0,1,2,... let h n := \p\ — 2/(tt\x\) + nxr c i n £ 2 (K 3 ), where Xr c 
denotes the characteristic function of the complement of T. Similarly, let : = 
\p + A\ — 2/(7r|x|) + nxr c - The diamagnetic inequality ()2.4|) and standard approx- 
imation arguments using Trotter's product formula imply that, for any u € L 2 (IR 3 ), 

|(exp(— th^)u)(x)\ < (exp(— th n )\u\)(x) . 

(See |FLSel Section 6.2] for details of the argument.) 

By the monotone convergence theorem |S1| Thm. 4.1], exp(— th n ) converges strongly 
to exp(— tH-p) on the subspace £ 2 (r), and similarly for h^. It follows that, for any 
ueL 2 (T), 

|(exp(-t^)M)(a;)| < (exp{-tH r )\u\){x) . 

Theorem 2.13 in [SHj yields the inequality || exp(— tH^ )|| 2 < || exp(— tH r )\\ 2 for the 
Hilbert-Schmidt norm, and hence || exp(— 2^^)11! < || exp(— 2tHr)\\i for the trace 
norm by the semigroup property. This holds for all t > 0, and hence proves (|4.3jl . ■ 

We use the notation (x)_ = max{0, —x} for the negative part of x e M. in the 
following. 

Lemma 4.2. Assume that there is constant M > such that 

Tr L2(r) (H r - A)_ < MA 4 (4.4) 

for all A > 0. Then 

Tr L2(r) (iJ r A -A)_<^MA 4 (4.5) 
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for all A > 0. 

We note the the numerical factor in (J4.5)) equals 6(e/4) 3 « 1.883. This factor is the 
price we have to pay, using our methods, to include an arbitrary magnetic field. It is 
the reason of the decrease of a c from 1/94 to 1/133. 

Proof. Since < e~ x ~ 1 ) we have 

Tr L2(r) (H£ - A)_ < ^Tr L2(r) exp {-tHf) 

for any t > 0. Using the diamagnetic inequality ()4.3j) . 

Tr L2(r) exp (-tHf) < Tr L 2 (r) exp (-tH r ) . 
Moreover, integrating by parts twice, e~ tx = t 2 J °° e _tA (x — A)_ dX, and hence 

POO 

Tr L2(r) exp (-t# r ) = t 2 / e~* A Tr L2(r) (# r - A)_ dX . 

Jo 

Using the assumption (|4.4jl . we thus obtain 

Tr L2(r) (H$ - A)_ < — — M / e- u A 4 dA = 24^M . 

e Jo t e 

To minimize the right side, the optimal choice of t is £ = 4/ A. This yields (|4.5jl . ■ 

In [LY, Thm. 11] it is shown that f)4.4j) holds for T = a ball of radius R centered 
at the origin. More precisely, the following proposition holds. 

Proposition 4.3. For any R > and A > 0, 

Tt lHBr) (H Br - A)_ < 4.4827i? 3 A 4 . 

Proposition 14.31 follows from Theorem 11 in |LY| by choosing \ to be the charac- 
teristic function of the ball Br, q — 1 and 7 to be the projection onto the negative 
spectral subspace of Hg R — A. 

Remark 4.4. It is illustrative to compare Proposition l4.3l with the Berezin-Li-Yau type 
bound 

TW, (Mr - A). < ^ I - A), to* = JLa'IFI . (4.6) 

(This can be proved in the same way as pVLol Thm. 12.3].) The right side of (jOj) 
is the semi-classical phase-space integral. The operator \p\ T is defined as H r above, 
but without the Hardy-term 2/(7r|x|). If the Hardy term were added, the phase-space 
integral would diverge (provided T contains the origin), but Proposition 14.31 savs that 
a bound of the form ()4.6)1 still holds. (An examination of the proof in |LY| shows that 
Proposition 14.31 actually holds for any open set T of finite measure.) 

Combining Lemma f4. 21 and Proposition 14.31 we obtain the following theorem, which 
replaces |LYt Thm. 11] in the magnetic case. 
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Theorem 4.5 (Lower bound on the short-range energy in a ball). Let C > 

and R > and let 

be defined on L 2 (R 3 ) as a quadratic form. Let < 7 < q be a density matrix (i.e., a 
positive trace-class operator) and let x by any bounded function with support in Br. 
Then 

Trxi X H£ R > -8.4411^|| X ||L • (4-7) 



As compared with |LY[ Thm. 11], the constant has been multiplied by 6(e/4) 3 , and 
llxllL appears instead of |-B_r| _1 ||xII2- 

Proof. Note that 

TrxixH$ R = Tr xiX - C/R) > -\\xix\U Ti lHBr) (h£ r - C/R). . 
The assertion follows from Lemma \4. 21 and Proposition ^..^ observing that ||x7X||oo < 

qWxWIo- ■ 

5. Proof of Theorem 12.11 

We assume that the reader is familiar with the proof of Theorem 2 in |LY| . We shall 
only emphasize changes in their argument. The main idea is to replace Theorems 10 
and 11 in |LY| by our Theorems 13.21 and 14.51 respectively. 

There are some immediate simplifications. First, in view of the simple inequality 
a/|p| 2 + m 2 > \p\ it is enough to prove the Theorem 12.11 for m — 0. Moreover, by 
the convexity argument of [DLJ it suffices to treat the case Z x = . . . = Z K =: z and 
az = 2/tt. So henceforth we assume m = 0, Z 1 = . . . = Z K = z and az = 2/ir. 

Let .Dfc := min{|i?^ — Ri\ : I ^ k} and define the Voronoi cell 

T fc := {x G M 3 : \x - R k \ < \x - R t \ for all I ^ k}. 
Fix < A < 1 and define a function W := G + F in each Voronoi cell by 

G(x) :=z\x-R k \-\ F(x):=D^F(\x-R k \/D k ), xeT k , 

where 

, 2- 1 (l-t 2 )- 1 if t < A 

b (t) : 



; v / 2i+|)r 1 if t > A. 

By the electrostatic inequality in |LY| Sect. Ill, Step A] our Theorem 12.11 will follow 
if we can prove that 

Tr j(\p + A\- aW) > £ (5.1) 

k=l 

for some < A < 1 and all density matrices 7 with < 7 < q). Note that ()5.1|) is an 
inequality for a one-particle operator. 
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For fixed < o < 1/3 we choose x, h as in (3.22), (3.24) in [LY]. Note that 
supp x C B\- a . Let 

Xk(x) := x{\ x ~ Rk\/Dk), h k (x) := h(\x - R k \/D k ). 
After scaling and translation, Proposition 13.21 yields that for any < 7 < q 
Tr 7 (|p + A| -aW) >TrxiJXi(\p + A\ -U^-aW) 

+ Tr(l - x?) 1/2 7(l - X\) 1/2 {\P + A\- W lie - aW) (5.2) 
- e^qn/Dx . 

Here, U^x) := D^U*{(x - i?i)/£>i) and Q, U* were defined in (J23J), (j3~B]l . (Note 
that our f2 is denoted by f2i in |LYj ) . Recall that [/* and Q are independent of A. 
We turn to the first term on the right side of ()5.2j) . Let C be a constant such that 

C>(l-a)(aF(\x\) + U*(x)^j for|x|<l-<7. (5.3) 

Note that xi is supported on a ball of radius (1 — a)Di centered at Ry. Hence 
aW(x) = (2/ix)\x — R\\~ l + D^F^x — Ri\/Di) on the support of Xi an d we can 
apply Theorem 14.51 to obtain the lower bound 

Tr X i7Xi(|£> -A\- Ul £ - aW) > Tr X ilXi (\p + A\ - 2 - C . n ) 

\ 7r|x — /ti| (1 — o)D\) 

aC A 

>-- M ^^h)Dr- (5 ' 4) 
We used also that |xi| < 1. Inserting (|5.4|) into (|5.2j) we find 

Tr 7 (|p + A| - aW) > -g.Dr 1 ^ 

+ Tr(l - x?) 1/2 7(l " X\) l/2 {\P + A\- Ul £ - aW) 

with 

A:= - + 8.4411- 



e (I -a) 

This estimate is exactly of the form (3.26) in |LYj . except for the value of the constant 
in A (which is called A in |LYj ). Starting from there one can continue along the lines 
of their proof. We need only note that in order to bound the last term in (3.29) in 
|LYj one uses the Daubechies inequality [Dj, which holds with the same constant in 
the presence of a magnetic field. (This is explained, for instance, in |LLoSi[ Sect. 5].) 
We conclude that stability holds as long as 



aq(A + J)<^, (5.5) 



where, as in |LY| Eq. (3.31)], 



J := 0.0258 



x\>l— 3(7 



2 

" +aF{\x\) + U*{x) 



7T\X\ 



dx . 
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This completes our proof of Theorem 12 .![ except for our bound on the critical a, which 
we justify now. 

As in |LYj . we choose a = 0.3, e = 0.2077 and A = 0.97. Our goal is to prove 
stability when qa < 1/66.5. We may assume a < 1/47, which is the assumption used 
in |LY| . Hence we can use the estimate J < 1.64 from |LY[ Eq. (3.40)]. 

To bound A, note that = 0.5571 as in |LY[ Eq. (3.30)]. It remains to choose 
an appropriate C satisfying ()5.3|) . For \x\ < 0.7 we have |a:| ) | < 1/1.02. Moreover, 
for U* we use the same estimate as in [LY], namely U*(x) < 0.2077 + 0.5751 = 0.7828. 
Using a < 1/(66.5 q) < 1/66.5, ()5.3|) therefore holds with 

0.7 (1/(66.5 ■ 1.02) + 0.7828) < 0.5583 =: C . 

This leads to a value of A = 1.7287. Hence ()5.5|) holds for qa < 1/66.5. 
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